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1. Introduction 

Since tlie discovery of tlie pure spinor formalism [Q the computation of manifestly 
supersymmetric superstring scattering amplitudes became possiblea. At first the results 
were limited to tree-level, where it was shown that amplitudes with an arbitrary number 
of bosonic and up to four fermionic states were equivalent to the standard results from 
the Ramond-Neveu-Schwarz (RNS) and Green-Schwarz (GS) formulations [^. Explicit 
four-point tree computations were first performed in 0, while in those results were 
streamlined in a superspace derivation which also made manifest its relation with one- 
and two-loop amplitudes. The five-point amplitude was computed in Q, providing a 
compact superspace representation which contrasts with the bosonic-only result from ||^. 
In addition, an OPE identity related to the Bern-Carrasco- Johansson kinematic relations 
I Tod was uncoveredp, which led to further developments discussed in [l^]. Furthermore, 



it was shown in |jT^ that there is a BRST-equivalent superspace expression for the field 
theory limit of the superstring amplitude from [^ which provides hints of a direct mapping 
between Feynman diagrams with cubic vertices and pure spinor superspace expressions. 

After being extended in [|T5| , p!6[| , higher- loop amplitude computations using the pure 
spinor formalism also became a reality. At one-loop, the massless four-point ||T5|jr^ , five- 



point [T^] and the gauge variation of the six-point amplitude [^ were obtained. At 



two- loops, the massless four-point amplitude was computed in [p0| , |21[] . Using formulae for 
integration over pure spinor space, the overall coefficients of the one- loop |^2| and two-loop 
p^ were also computed and shown to agree with S-duality conjecture expectations [|2^ . 



Besides streamlining amplitude computations avoiding supermoduli spaces and sums 
over spins structures altogether, the pure spinor formalism naturally gives rise to manifestly 
supersymmetric kinematic factors in pure spinor superspace. 

Pure spinor superspace expressions are correlation functions written in terms of ten- 
dimensional super- Yang-Mills superfields [^| and three pure spinors A" normalized agj 



{{x-f'^e){Xj^e){XYe){0jmnp0)) = i. (i.i) 

The simplest example of a pure spinor superspace expression is provided by the massless 
three-point scattering amplitude [|l[ 

irf = {{XA'){XA'){XA^)). (1.2) 



For reviews see 

See also [11,12[ for string theory monodromy explanations of the BCJ relations. 



The precise overaU coefficients of [2^,23] wiU not be needed here. 



The four-point kinematic factors are given by 



i^f = {{\A^){x^^w''){\rw^):Ftr:), 



K 



(2) 



{{Xi^'^^'^^x){\Yw^):fI^:fI:fI) 



;i.3) 

;i-4) 
:i.5) 



mn'^ pq"^ rsi 

for the tree-level [0 , one- |]15| and two-loop ( |1.5| ) amplitudes, respectively. Another example 
is provided by the one- loop five-point kinematic factor of [|18[, whose expression for the 
(12) and (25) "channels" read 



Li, 



{XA^){k^ . A^) + AMXYW^) {X^^W'){X^-W')Ft^) 
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K25 = ((AAi) {x^^w^){e ■ A'^) - -{xty'w'>):fI (A7"T^^)J-;^J-(2^5). 



The above pure spinor superspace representations provide compact information about 
the amplitudes, but it may be convenient to evaluate those expressions in terms of familiar 
component expansions. These component expansions are written in terms of polarization 
vectors e^ and spinors xj with momenta /c^, where / = 1, . . ., A^ are the particle labels and 
m = 0, . . ., 9 a = 1, . . ., 16 are the Lorentz and Weyl indices of ten dimensional Minkowski 
space. 

The general method to evaluate these expressions in components was explained in the 
appendix of [^ . One uses the 6'-expansions |^,|| of the S YM superfields 



o i/ 
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;i.6) 



where am{x) = e^e*'^^ and ^"(x) = x'^e*'^^. After that, only terms containing five ^'s are 
kept. Using symmetry alone it is possible to rewrite arbitrary combinations of {X^9^) in 



terms of Kronecker deltas and epsilon tensors [pT| , p7| . For example 

1 



((A7"^^)(A7"^)(A7^^)(^7..s^)) 

2 



Xmnp 



;i.7) 



Substituting the various correlators by their corresponding tensors as above, the component 
expansion of pure spinor superspace expressions is obtained. 

Several different computer-aided procedures were used along the past years for the 
above steps in superspace derivations of scattering amplitudes. At first, those pure spinor 
correlators of W^ were obtained with the helpB of the GAMMA package [EH]. With some effort. 



the superspace expressions were expanded in ^'s by hand, the corresponding correlators 
identified from a catalog composed of entries like (|1.7|) and the resulting tensors were typed 
in Mathematica or FORM |j2^ . When the scattering involved fermionic states, various Fierz 
identities were usually necessary at intermediate stages. After all these steps, the output 
consisted of several terms composed of Kronecker deltas and Levi-Civita epsilon tensors 
contracted with momenta and polarizations. Those terms constituted the final answer. 

However the semi-automated method described above does not scale very well for 
higher-point amplitudes or when there are many expressions to evaluate in sequence. This 
fact provided the motivation to write the program presented here. It was developed mainly 
to help the author's own workfiow during computations, and therefore it reflects his pri- 
orities (and it is a continuous work in progress). It is called PSS and it is written in the 



interesting language of FORM [p9 |. 



2. How PSS works 

The goal is to be able to obtain component expansions in a fully automated process - 
all that is required is the pure spinor superspace expression to be expanded and whether 
the external states are bosonic or fermionic. The rest must be done by the computer. 



For example, the superspace expression for [^ is typed in PSS as 



Local [4-pts_one-loop] = la*Al*la*ga(in)*W2*la*ga(n)*W3*cF4(m,n) ; 



A small C program was also written to deal with anti-symmetrization of huge tensors during 
consistency checks. 



Choosing all fields to be bosonic results in the following ipsis litteris output, 

[4-pts_one-loop] = 

- l/5760*kl.e2*kl.e3*el.e4*t 
+ l/5760*kl.e2*kl.e4*el.e3*t 
+ l/5760*kl.e2*k2.e3*el.e4*u 
+ l/5760*kl.e2*k2.e4*el.e3*t 
+ l/5760*kl.e2*k3.el*e3.e4*t 
+ l/5760*kl.e3*k2.el*e2.e4*t 

- l/5760*kl.e3*k2.e4*el.e2*t 

- l/5760*kl.e3*k3.e2*el.e4*u 

- l/5760*kl.e3*k3.e2*el.e4*t 
+ l/5760*kl.e4*k2.el*e2.e3*u 

- l/5760*kl.e4*k2.e3*el.e2*u 
+ l/5760*kl.e4*k3.e2*el.e3*u 
+ l/5760*kl.e4*k3.e2*el.e3*t 

- l/5760*k2.el*k2.e3*e2.e4*u 
+ l/5760*k2.el*k2.e4*e2.e3*u 
+ l/5760*k2.el*k3.e2*e3.e4*u 

- l/5760*k2.e3*k3.el*e2.e4*u 

- l/5760*k2.e3*k3.el*e2.e4*t 
+ l/5760*k2.e4*k3.el*e2.e3*u 
+ l/5760*k2.e4*k3.el*e2.e3*t 
+ l/5760*k3.el*k3.e2*e3.e4*u 
+ l/5760*k3.el*k3.e2*e3.e4*t 
+ 1/I1520*el.e2*e3.e4*t*u 

- 1/I1520*el.e3*e2.e4*t*u 

- 1/I1520*el.e3*e2.e4*t"2 

- 1/I1520*el.e4*e2.e3*u"2 

- 1/I1520*el.e4*e2.e3*t*u 

i 

Momentum conservation: k4 eliminated 
Gauge invariance: not tested 

0.05 sec +0.08 sec: 0.13 sec out of 0.15 sec 



which is the result obtained in [O]- One should notice in the final statistics displayed by 



FORM how quickly the answer is obtained. 

The program is composed of one main FORM script called pss .frm and four header 
files: pss_header .h, kin_f actor .h, pss.h and ps_tensors.h. They contain the defi- 
nitions of indices, vectors, tensors, superfields etc (pssjieader .h), the pure spinor su- 
perspace expressions to be evaluated (kin_f actor .h) and the procedures which actually 
do the computations (pss.h). The database of pure spinor correlators is contained in 
ps_tensors.h. There is also a small sed script (F0RM2tex.sed) to help translating the 
result into T^. 

To use the program one has to write down the kinematic factor in kin_f actor .h. In 
the beginning of the main file pss . frm, the number of points must be defined (e.g. #def ine 
Npts ' '4' ') and whether the external states are bosonic or fermionic (#def ine f ieldl 
' ' ' ' and so forth). After that one executes pss . frm using either form or tform (for multi- 
processor computers) and a result like the one written above is obtained. Optionally one 
can select which momentum to be eliminated by setting the dollar variable $kn and test 
for gauge invariance by uncommenting the line containing id el = kl; #gauge = el; 
close to the end of the file. Several other things can be done, depending on the problem 
at hand and how one chooses to manipulate it. There are also a few debug options (-d 
psonly, sfexpand, nofierz) which help in case something goes wrong and one has to check 
where. They will be explained below. 

2.1. User input and notation 

The superspace expression is written in terms of the super- Yang-Mills superfields 
A^(6', x), Al^{6, x), W^"(6', x) and J^^ni^^ ^) l^- Their definitions are contained in the file 
pss_header.h and they correspond to, 

Al, Bl, Wl, cFl(m,n) 

The pure spinor A" and 6" are denoted by la and th. Note that PSS does not know 
about spinor indices, but that does not cause problems as long as one writes down correct 
superspace expressions. However, that also means one has to take care whether a fermionic 
superfield is contracted from the left or right, because the ^-expansions to be used differ 



in this caseH. For this situation one has to use the "left" version of the superfields, ALl 
and WLl etc. For example, the factor (W^'jrnnpW^) must be written as 

WLl * ga(m, n, p) * W2 

The 16 X 16 gamma matrices 7"^, 7"^", j-mnp^ ^^^^ denoted by 

ga(m), ga(m,n), ga(m,n,p), ... 

Another convention to be followed is to write the pure spinor A" to the left of fermionic 
bilinears, so {X^'^W'^) or {W'^j'^X) must be written as 

la * ga(m) * W2 

and not 

WL2 * ga(m) * la 

although PSS can be easily modified to accept the latter version. Of course, there is no 
problem to write the factor (X^'^^^'^^X) and the procedure to identify correlators is aware 
of it. The generalized Kronecker delta is defined byl3 

^'C^-.nT = da(mi, . . ., mAT, ni, . . ., n^) 

2.2. The computation 

It is important to understand how PSS actually obtains the component expansions, 
so that modifications can be easily done. All the action happens inside pss . f rm, where it 
calls the procedures from pss.h. Let us now follow some of the stepso. 

When PSS is executed it loads the headers and the kinematic factors. The kinematic 
factors are local variables which will be manipulated by the FORM program. The first 
part of the manipulations transform the superspace input into an expression suitable for 



This is an artifact of how PSS was designed and this distinction is meaningless in real life. 
Unfortunately FORM has no notion of a generahzed Kronecker delta, so it had to be defined 
by PSS. Note that the usual Kronecker delta is defined by FORM as 5^ = d_(m,n). 
One should also read the source code, as some details will be skipped. 
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identifying the required pure spinor correlators in terms of tensors to arrive at the final 
answer. The (slightly simplified) sequence of procedures is the following: 

#call chooseMomentumO 

#call superf ieldExpandO 

#call keep5thetas() 

#call gammaExpandO 

#call deltaExpandO 

#call onShellO 

#call orderFermionsO 

#call PSorderingO 

#call FierzO 

#call earlySimplifyO 

#call identifyCorrelatorsO 

#call deltaExpandO 

#call dualizeGammasO 

#call fieldStrengthO 

#call onShellO 

#call gammaExpandO 

#call orderFermionsO 

#call diracEquationO 

#call momentumConservationO 

These steps are almost self-explanatory, and they correspond to what one would actually 
do in a computation with pen and paper. 

The procedure chooseMomentumO is called to choose which momentum to eliminate. 
Although the result does not depend on the choice, there are differences in the number 
of intermediate terms and computing time. For example, when there are fermionic fields 
it is convenient to eliminate one momentum whose label is not one of the labels of the 
fermionic particles. The reason is the increased chance of applying the Dirac equation on 
fermion bilinears to reduce the rank of the gamma matrix. For example, {x^l^^^')(?)^m — 

2(xW[V)fc?- 

Then the procedure superf ieldExpandO is called, which expands the SYM super- 
fields in their 6 components using ( |1.6| ) and selects the terms according to whether the 



particles are bosonic or fermionic. The assignment of particles to each label of the super- 
fields is done at the beginning of pss .frm with 

#define fieldl "0'' 
#define field2 "0'' 



where (1) means bosonic (fermionic). Terms which do not contain five ^s are then 
discarded with keep5thetas(). At this point, one of the terms in ( |1.4| ) looks like 

- l/256*la*ga(el) *th*la*ga(k4) *ga(m2 ,n2) *th*la*ga(n) *ga(m3 ,n3) *th* 
th*ga(n,m4,n4)*th*F2(m2,n2)*F3(m3,n3)*F4(m4,n4) 

which is the FORM output for 



256 

in the schoonschip notation. The gamma matrices in (p.l| ) are expanded using gammaEx- 
pandO, 

^n mana _ nmzn^ _|_ man 713 _ n^n m^ 

More complicated expansions may introduce generalized Kronecker deltas, which are then 
expanded with deltaExpandO, e.g., S'^J' = ^{6^6^ - S^6^). The procedure onShellO 
kills any term which may have been generated at this point containing {k^ ■ k^) or (/c* • e*). 
When there are fermionic external states, the procedure orderFermionsO rewrites 
the fermionic bilinears (;y*7"^^'"^"X-') such that i < j. This procedure keeps track of 
overall minus signs which may be needed due to the Grassmanian nature of x's and the 
symmetry properties of the gamma matrices. For example [x^l^^X^) ~^ ~(x^7"^X"^) o^^ 

(x'7"'"V)^+(xW™V)- 

The procedure PSorderingO follows a set of conventions on the ordering of fermionic 
bilinears to minimize the number of pattern matching when trying to identify the pure 
spinor correlators needed for one particular computation. For example, if the expression 
{\'^'^6){6'^qrsd){,^l'^d){,^l^d) is encountered in the middle of the computations, it is first 
rewritten by PSorderingO as {X'j^9){Xj'^9){X'yP9){9'yqrsd) ■ Later on there will be only 
one pattern to match in order to identify the correlator (in identifyCorrelatorsO) and 
replace it with ^da(m,n,p,q,r,s). For the same reason, the ordering "inside" gamma 



matrix bilinears is such that the fields appear as (A7"^'''6'), (A7"^'x) and (x7'^"'6'), and 
not for example as (^7'^'A) etc. 

When the computation involves fermionic particles, after processing the expressions 
with PSorderingO there may be factors such as (x^7m^)(x^7'^^), which will be rewritten 
by the procedure FierzO as 

-^(xS"^7.r.7"x')(^7^^^^) 

again in order to minimize the number of pattern matching to identify correlators. This 
procedure takes care of various bilinears combinations, and the pair of fermions chosen to 
be expanded depend on the particular combination being considered, for which the user 
should read the source code for more information. The general expansions are done with 
the formula 

A-x'^ = ^(A7^X)7;^'' + ^i^l'^-'^'xh^, + ^{^'^''''''^Xh^p.r- (2.2) 

If the procedure FierzO takes effect, then there will be more gamma matrices to expand, 
like in the example above (x^7"^7grs7"^X^)- If that is the case, then the extra calling of 
gammaExpandO and deltaExpandO, followed by PSorderingO and orderFermionsO 
will let the expressions ready for being identified. 

Using the example of the one-loop kinematic factor, at this point one of the terms 
being dealt with by PSS is 

- l/512*la*ga(Nl_? ,N2_? , N3_?) *th*la*ga(N4_? , N5_? ,N6_?) *th* 
la*ga (N7_?) *th*th*ga (N7_? , N8_? , N9_?) *th*Fl (N8_? , N9_?) * 
F2 (Nl_? , N2_?) *F3 (N4_? , N5_?) *F4 (N3_? , N6_?) 



which is 

1 



512 



(A7-^'^^-^^)(A7™"«^)(A7-^^)(07"^-«"^^)F4„^F^2^^^F4„^F4„^. (2.3) 



The next step is to call the procedure which identifies the pure spinor correlator from a 
catalog of known tensors. This is done with identifyCorrelatorsO, after which the 
above term is given by 

- 1/512*F1(N1_?,N2_?)*F2(N3_?,N4_?)*F3(N5_?,N6_?)*F4(N7_?,N8_?)* 
ps331(N3_?,N4_?,N7_?,N5_?,N6_?,N8_?,N9_?,N9_?,Nl_?,N2_?) 

9 



That is, PSS identifies the pure spinor correlator {X-f''^''^''3 0^(^x^nin5ne0^(^^^m0^(^0^nrnsnQ0^ 
with the tensor ps331(?2i, n2, ns, n^, n^, ug, nj, n^, ng, ng). Looking at the appendix of |1T9[] , 
this tensor is expanded in terms of Kronecker deltas as 

{{X^^^Pe){Xj'^^'9){X-ft9)i9-f,jkO)) = -l-e^^'=™P9-*+ (2.4) 



1 



140 



"t "[i'l '^j'^k] "t 0[i'l "j "fcjj 280 I ' -^ 1 ^ -^ ^^ 

and this is one of the correlators included in the catalog ps_tensors .h. So the purpose of 
identifyCorrelatorsO is to transform an input containing correlators with {X^9^) into 
an expression written in terms of tensors like ps331(m,n,p,q,r,s,t,u,v,x), which are later 
substituted by Kronecker deltas and epsilon tensors as in ( |2.4|) using the catalog. 
The list of correlators identified by identifyCorrelatorsO is 

((Af^-*^A)(A7'"^)(^7'^-'^^^)(^7^^-^^^))=uind(ti,...,t5,m,ri,...,r3,.i,...,S3) 

((A7*-*^A)(A7"^-"^^^)(^7'^-'^^^)(^7^-^^^)) = tind(ti,...,t5,mi,...,m3,ri,...,r3,si,...,S3) 
((A7*-*sA)(A7"^i-"^^^)(^7""-'^^^)(^7'"-"^^)) = cind(ti,...,t5,mi,...,m5,ri,...,r3,si,...,S3) 
{{X'y'^9){X'y''9){XY9){9-f''^''9)) = pslll(m, n,p, a, b, c) 
{{X-f'^9){X-f''9){X'yP9){9'y'''-'^'9)) = psllleps(m,7i,p,ri, ...,r7) 
((A7*-*^^)(A7"^^)(A7"^)(^7'^^-'^^^)) = ps311(ti,...,t3,m,n,ri,...,r3) 
((A7*^-*^^)(A7"^^)(A7"0)(^7'^i-'-^^)) = ps311eps(ti,...,t3,m,n,ri,...,r7) 
((A7"^-"^^^)(A7--"^^)(A7"^)(^7""-'^^^))=ps331(mi,...,m3,ni,...,n3,a,ri,...,r3) 
((A7"^i-"^^^)(A7"i-'^^^)(A7"e)(e7""-"^^))=ps331eps(mi,...,m3,ni,...,n3,a,ri,...,r7) 
((A7*^-*3^)(A7"^-"^^^)(A7"-"^^)(^7'^i-^^^))=ps333(ti,.,t3,mi,.,m3,ni,.,n3,ri,.,r3) 
((A7*-*30)(A7"^-"^^^)(A7--"3^)(^7'^-'^-^))=ps333eps(ti,.,t3,mi,.,m3,ni,.,n3,ri,.,r7) 
((A7*i-*^^)(A7"^e)(A7"^)(^7""-"'^))=ps511(ti,.,t5,m,a,ri,.,r3) 
((A7*^-*^e)(A7"^^)(A7"^)(^7^-'^^^)) = ps511eps(ti, ., t,, m, a, n, ., rr) 
((Afi-*^^)(A7"^-"^3^)(A7"^)(^7^-'^^^))=ps531(ti,.,t5,mi,.,m3,a,ri,.,r3) 
((A7*^-*^^)(A7"^i-"^^^)(A7"^)(^7''^-'^^^))=ps531eps(ti,...,t5,mi,...,m3,a,ri,...,r7) 
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((A7*^-*^0)(A7'--"^^^)(A7"--30)(^^.i...r30))=pg553(^^^.^^^^^^^.^^^^„^^.^^3^^^^.^^3) 

((A7*^-*^^)(A7"^-"^^^)(A7"-"^0)(^7'^-'-^0))=ps533(ti,.,t5,mi,.,m3,ni,.,n3,ri,.,r3) 

It may happen that one particular computation requires a correlator not in the list. When 
that happens PSS automatically detects the missing correlator and prints it before exiting. 
For example 

I do not identify a correlator for: 
+ l/983040*chil*ga(Ml_?,N2_?,N3_?,N4_?,N5_?,N6_?,M7_?)*chi2* 
la*ga (Nl_? , M2_? , N3_? , N4_? , N5_?) *th*la*ga (kl , N6_? , N8_?) *th* 
la*ga(N9_? , N10_? , Nll_?) *th* 

th*ga(N7_? , N8_? , N9_? , N12_? , N13_? , N14_? , N15_?) *th* 
F3(N10_?,N11_?)*F4(N12_?,N13_?)*F5(N14_?,N15_?) 
Add it to the identifyCorrelatorsO procedure 

In this case the missing correlator can be obtained from ps533() and its pattern matching 
added to the procedure identifyCorrelatorsO and its tensor representation added to 
the file ps_tensors.h. 

After calling identifyCorrelatorsO the tensors above are substituted by their ten- 
sor representations. This is done by including the file ps_tensors .h, which contains their 
expansions in terms of generalized Kronecker deltas and Levi-Civita epsilons. The proce- 
dure deltaExpandO expands the generalized Kronecker deltas in terms of the antisym- 
metric combinations of the usual Kronecker delta. Gamma matrices with more than five 
indices are manipulated with dualizeGammasO, where for example the following identity 
is used 

mi. ..my _ __„v ^ni...n3 

lap „| '-Cmi...m7ni...n3 fa/? 

The cases where epsilon tensors are contracted with gamma matrices is also dealt with, 
For completeness, the gamma matrix conventions [^ are such thatB 

mi...mg ■ mi...mQn 



^■■-i o _ , p..«i -y-if'^ "l „ (^rni...ms\ P _ • mi...m8ni ...n2 / \ P 



2! 



The signs change when both spinor indices of the matrix matrices change from Weyl to 
anti-Weyl. 
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mi. ..my ^- m-i_...mrni...na ni...n:i ( mi...me\ p _ _}_ • mi...m6ni...n4 / N /3 



31"- icxp ' V/ /a ^1 

■ mi...m5ni...n5/ N ^mi...m3__J_ 



Tni...Tn5 _ J_ ■ mi...m5ni...n5/ N ^mi...m3 _ J_ ■ mi ...m3ni...n7 / \ 

'a/3 ~ K} Wni...n5ja/3, 7a/3 — y| * t (,7".! ...^7 Ja/3 



(2.5) 

It is important to notice that FORM uses the convention that emi...mio^^^"'^^° = 10! instead 
of —10!, so that is why there are factors of i together with epsilon tensors in PSS. 

The procedure f ieldStrengthO substitutes F^^ = fc^e^ — k^^^^ ^ind onShellO 
annihilates terms with (k^ ■ k^) and {k^ • e^). If there are fermionic particles, the procedure 
diracEquationO uses the Dirac equation to reduce the rank of the gamma matrices when 
there is a momentum of one of the particles in the fermionic bilinear being contracted with 
one of its indices, for example 

(xWmnfcaX^) = A;^(xWnX^) " ^^(xWmX^)- 

Finally, the procedure momentumConservationO applies the conservation of momentum 
to one of the labels, which can be manually chosen by setting the "dollar" variable $kn in 
the beginning of pss.frm (if let at its default value of zero, then an automatic choice is 
made) . 

After the above (simplified) sequence of steps the desired component expansion of the 
superspace expression is printed on the screen. 

If one chooses the particles 1 and 2 to be fermionic by using #def ine f ieldl 1 etc, 
rerunning the program results in, 

[mafra(9Pilar:pss] tform -q -w2 pss.frm 
[4pts_one-loop] = 

+ 1/I1520*chil*ga(k4, e3 , e4) *chi2*u 
+ 1/I1520*chil*ga(k4, e3 , e4) *chi2*t 

- 1/I1520*chil*ga(k4)*chi2*e3.e4*u 
+ 1/I1520*chil*ga(k4)*chi2*e3.e4*t 
+ l/5760*chil*ga(e3)*chi2*kl.e4*u 

- l/5760*chil*ga(e3)*chi2*k2.e4*t 
+ 1/I1520*chil*ga(e4)*chi2*k4.e3*u 

- 1/I1520*chil*ga(e4)*chi2*k4.e3*t 

Momentum conservation: k3 eliminated 
Gauge invariance: not tested 

0.08 sec +0.42 sec: 0.50 sec out of 0.32 sec 
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2.3. Debug options 

There are three pre-defined debug options which can be used to check intermediate 
steps in the computation, sf expand, psonly and nof ierz. If sf expand is invoked (using 
the -d fiag; see FORM's manual), then only the superfield expansion in terms of ^'s is printed, 

[mafra(9Pilar:pss] tform -q -w2 -d sf expand pss.frm 

[4pts_one-loop] = 

- l/256*la*ga(el) *th*la*ga(k4) *ga(m2 ,n2) *th*la*ga(n) *ga(m3 ,n3) *th* 
th*ga(n ,m4 , n4) *th*F2 (m2 , n2) *F3 (mS , n3) *F4 (m4 , n4) 

- l/384*la*ga(el) *th*la*ga(m) *ga(k2 ,N1_?) *th*th*ga(m2 ,n2 , Nl_?) *th* 
la*ga(n) *ga(m3 , n3) *th*F2 (m2 , n2) *F3 (m3 , n3) *F4 (m , n) 

+ l/256*la*ga(el) *th*la*ga(m) *ga(m2 ,n2) *th*la*ga(k4) *ga(m3 ,n3) *th* 

th*ga(m ,m4 , n4) *th*F2 (m2 , n2) *F3 (m3 , n3) *F4 (m4 , n4) 
+ l/384*la*ga(el) *th*la*ga(m) *ga(m2 ,n2) *th*la*ga(n) *ga(kl , Nl_?) *th* 

th*ga(m3 , n3 , Nl_?) *th*F2 (m2 , n2) *F3 (m3 , n3) *F4 (m , n) 
+ l/384*la*ga(el) *th*la*ga(m) *ga(m2 ,n2) *th*la*ga(n) *ga(k2 , Nl_?) *th* 

th*ga(m3 , n3 , Nl_?) *th*F2 (m2 , n2) *F3 (m3 , n3) *F4 (m , n) 
+ l/384*la*ga(el) *th*la*ga(m) *ga(m2 ,n2) *th*la*ga(n) *ga(k4,Nl_?) *th* 

th*ga(m3 , n3 , Nl_?) *th*F2 (m2 , n2) *F3 (m3 , n3) *F4 (m , n) 

- l/512*la*ga(Nl_?)*th*th*ga(ml,nl,Nl_?)*th*la*ga(m)*ga(m2,n2)*th* 
la*ga(n)*ga(m3,n3)*th*Fl(ml,nl)*F2(m2,n2)*F3(m3,n3)*F4(m,n) 

which can be useful to check when something goes wrong. The debug option psonly prints 
the expression after the correlators were identified, for example, 
[mafra(9Pilar:pss] tform -q -w2 -d psonly pss.frm 

[4pts_one-loop] = 

- 1/128*F1(N1_?,N2_?)*F2(N3_?,N4_?)*F3(N5_?,N6_?)*F4(N3_?,N5_?)* 
pslll(N7_?,N4_?,N6_?,N7_?,Nl_?,N2_?) 

- 1/256*F1(N1_?,N2_?)*F2(N3_?,N4_?)*F3(N5_?,N6_?)*F4(N3_?,N7_?)* 
ps311(N5_?,N6_?,N7_?,N8_?,N4_?,N8_?,Nl_?,N2_?) 

- 1/256*F1 (Nl_? , N2_?) *F2 (N3_? , N4_?) *F3 (N5_? , N6_?) *F4 (N5_? , N7_?) * 
ps311(N3_?,N4_?,N7_?,N8_?,N6_?,N8_?,Nl_?,N2_?) 

- 1/512*F1(N1_?,N2_?)*F2(N3_?,N4_?)*F3(N5_?,N6_?)*F4(N7_?,N8_?)* 
ps331(N3_?,N4_?,N7_?,N5_?,N6_?,N8_?,N9_?,N9_?,Nl_?,N2_?) 
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where the other terms are similar and were omitted. Furthermore, nof ierz prints the 
superspace expansions before any Fierz manipulation is done, 

[mafra@Pilar :pss] tform -q -w2 -d nof ierz pss.frm 
[4pts_one-loop] = 
+ l/24*cliil*ga(Nl_?)*th*chi2*ga(N2_?)*th*la*ga(k2,N2_?,N3_?)*th* 

la*ga(Nl_?) *th*la*ga(N4_?) *th*F3 (N4_? , N5_?) *F4 (N3_? , N5_?) 
+ l/48*chil*ga(Nl_?) *th*chi2*ga(N2_?) *th*la*ga(k2 , N2_? ,N3_?) *th* 

la*ga(N4_? , N5_? , N6_?) *th*la*ga(Nl_?) *th*F3 (N5_? , N6_?) *F4 (N3_? , N4_?) 
+ l/24*chil*ga(Nl_?) *th*chi2*ga(N2_?) *th*la*ga(Nl_?) *th*la*ga(k2) *th* 
la*ga(N3_?) *th*F3 (N3_? , N4_?) *F4 (N2_? , N4_?) 



There are many possible extensions and optimizations which can be made to PSS, as 
it is available to download at [http : //www . aei . mpg . de/~crmaf ra/pss . tar . gz| under the 
GPL license. In particular, dealing with four-fermion expansions is still not completely au- 
tomated (nor guaranteed to be correct) . It would be interesting to implement the fermionic 



methods described in \^7j for this purpose. Furthermore, it should be straightforward to 
write procedures to translate the full ten-dimensional components to four dimensions using 
the spinor helicity formalism, in order to compare with the results appearing in [^. The 
possibilities are many and it is hoped that PSS provides a framework for further work. 
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